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Within the framework of the (3+l)-dimensional Lorentz-violating extended electrodynamics in- 
cluding the CPT-odd Chern-Simons term, we consider the electromagnetic field between the two 
parallel perfectly conducting plates. We find the one-particle eigenstates of such a field, as well as 
the implicit expression for the photon energy spectrum. We also show that the tachyon-induced vac- 
uum instability vanishes when the separation between the plates is sufficiently small though finite. 
In order to find the leading Chern-Simons correction to the vacuum energy, we renormalize and 
evaluate the sum over all one-particle eigenstate energies using the two different methods, the zeta 
function technique and the transformation of the discrete sum into a complex plane integral via the 
residue theorem. The resulting correction to the Casimir force, which is attractive and quadratic in 
the Chern-Simons term, disagrees with the one obtained in [M. Frank and I. Turan, Phys. Rev. D 
^ , 74, 033016 (2006)], using the misinterpreted equations of motion. Compared to the experimental 

' data, our result places a constraint on the absolute value of the Chern-Simons term. 
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I. INTRODUCTION 

For today, the Standard Model is proved with a convincing variety of experiments. However, this theory cannot 
_ . shed light on some essential aspects of our understanding of the Nature, which still remain obscure and are of great 
interest. Among them is the quantum gravity, which is believed to play a major role at the Planck scale of energies 
(— I ■ {Epi ~ lO^^GeV). The Standard Model fails to incorporate the General Relativity at the quantum level, since the 
'— ' I quantization methods adopted in the former theory lead to a nonrenormalizable quantum gravity in this case. Such 
theory is unacceptable as a fundamental theory. At the same time, there exist some candidates for such Fundamental 
theory, string theories, for instance, taking the form of the Standard Model and the General Relativity in appropriate 
(*~^ low-energy limits. Thus searching for signatures of the Planck-scale physics at experimentally attainable energies 
QQ [ would be a natural way to choose between these candidate theories, to constrain their parameters, and to clarify the 

• essence of the quantum gravity. 
CO ' Planck energies being far from experimental attainment, the Standard Model Extension (SME) was elaborated. 

■ It is an effective theory (applicable at the energies E ^ E-pi) formulated axiomatically as a set of corrections to 

■ the Lagrangian of the Standard Model, which fulfill some "natural" requirements P, such as observer Lorentz 
0^ invariance, ^-inomentum conservation, unitarity, and microcausality. In what follows, we will focus on a subset of 

the SME referred to as the minimal SME in flat Minkowsky spacetime, in which local SU{^)c x SU{2)i x U{1)y 
^ \ gauge invariance and power-counting renormalizability are also required. A spectacular feature of such requirements 
is that they reduce the diversity of possible correction terms down to a finite number of them. Each of them consists 
of a complex (pseudo)tensor constant (SME coupling) contracted with conventional Standard Model fields and their 
i spacetime derivatives. These constants are believed to stand for vacuum expectation values of the fields featuring in 
■ " ' the hypothetic Lorentz-covariant Fundamental theory and condensed at low energies due to the spontaneous symmetry 
breaking mechanism. Indeed, it has been shown that such Lorentz symmetry breaking can occur in some theories 
beyond the Standard Model 0, H i, H, leading subsequently to the SME. The SME can thus be used to reduce 
the complexity of these theories and related calculations in the low-energy limit. It also provides a standard for 
representation of the data obtained in experiments searching for Lorentz violation. 

Recently, a number of theoretical researches has been performed aimed at investigating the vacuum structure of 
this model (see, e^., 0, S [13, Uli), and to study the assumed Lorentz violation on various high-energy processes 
[T^ [Tsl . [3, [3, [3- Such effects as the vacuum photon splitting, vacuum birefringence, and vacuum Cherenkov 
radiation were studied. The effect of Lorentz violation on the synchrotron radiation [itI [TsI and the atomic spectrum 
and radiation properties [l^, [l^, [2l| were also analyzed. As a result of such theoretical investigations and subsequent 
precise experimental tests, the new constraints on Lorentz violation were placed. 
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On the other hand, even in the conventional Standard Model, the quantum structure of the vacuum manifests itself 
in such effects as the Casimir effect [13, [23], that have been directly observed. This effect has been studied thoroughly 
within the Standard Model, including different approaches [13, [1^ (various regularization schemes, the calculations 
via the dyadic Green function and the vacuum energy), setups (two parallel plates, a sphere, a plate and a sphere; 
non-perfect conductor plates, finite temperature, fermion Casimir effect, etc.). The so-called Maxwell-Chern-Simons 
(MCS) electrodynamics in (2 + 1) dimensions was also considered in this concern [2^. The Casimir effect was also 
studied within the context of extra dimensions (see, e.g., [23,[H, [s^]) and curved space [1^, [3l|, [s^l . Some attention 
was paid to this problem within the SME, namely, within Extended Quantum Electrodynamics (see Sec. |TT]) [ssj . 
which, in a certain particular case, takes the form quite analogous to the MCS electrodynamics. 

Although the original version of Casimir effect concerned electromagnetic field, this effect is much more general. It 
can be defined as the stress (force per unit area) on the bounding surfaces when a quantum field is confined within a 
finite volume of space. The boundaries can be material media and also the interfaces between different phases of the 
vacuum or topologies of space. 

In our investigation, we focus on the extended electrodynamics, in order to find the signature of the Lorentz violation 
in the electromagnetic Casimir effect. 

Our paper is organized as follows. In Sec. |TT1 we analyze the pure-photon sector of a particular case of the SME, 
namely, of the (3 -I- l)-dimensional Maxwell-Chern-Simons electrodynamics, in order to obtain the expression for the 
vacuum energy within this theory. We prove directly that the Casimir force in such a theory remains gauge-invariant 
although the energy-momentum tensor does not. In Sec. IIIIl we focus on finding the eigenstates and the energy 
spectrum of the photon between the two parallel perfectly conducting plates within the MCS electrodynamics. The 
one-photon energy spectrum, which follows from the implicit expression (I3.22p found in Sec. Illli is then analyzed in 
Sec. IIVI where we show that when the plates are close enough to each other, the imaginary-energy (tachyonic) states 
are negligibly few, so the instability of the theory is also negligible (see assertion IH]) . In Sec. |Vl using approximate 
energy eigenvalues, we find the leading correction to the zeta function corresponding to the one-photon energy-squared 
operator, and, in turn, the leading correction to the Casimir energy. The contribution of the "quasi-zero" modes, 
which become trivial in the Maxwell case, to the Casimir force, is also discussed. In Sec. I VII we use the strict 
approach based on the residue theorem |34| to explicitly sum the vacuum energy series, which contains the zeros 
of the transcendental function (|3.22[) . We then renormalize the resulting complex plane integral and find the exact 
convergent expression (|6.20p for the real part of the Casimir energy. Section IVIII summarizes the results obtained, 
and we discuss the constraints they place on the parameters of the Lorentz violation. 



II. EXTENDED MAXWELL-CHERN-SIMONS ELECTRODYNAMICS 



A. General notes 



After the electroweak symmetry is broken within the minimal SME, the extended Quantum Electrodynamics arises 
[l[. For one fermion generation, it contains 10 additional (tensor) coupling constants which introduce the interaction 
between the Dirac electron, the photons, and the condensates of Planck-scale fields whose nonzero values violate 
Lorentz invariance at low energies: 

^ext.QED = -\^F^,F^''' + 4>{iC''D,- M)^~^{kF)^,,o.pF^'''F''f^ + r^^A''F^,, (2.1) 

r- = ^- + c^-^^ + dM-^5^^ + e- + ,/-^^^+ 1^«/3V,^^ (2.2) 

M = m + a,,7^ + 6^757^ + (2.3) 

where _F),^ = ^^/ivapF"^ is the dual field strength tensor, e^j/Q/j is the Levi-Civita symbol, with e^^^^ = — £0123 = +1, 
F>i, ~ d,y+ieAi, is the U (l)-covariant derivative, and m and e are the electron mass and charge, respectively. 7^, 75, and 
cr^y denote the Dirac matrices, while spacetime indices fi, i>,a, (3 = 0, 1, 2, 3. For the Lorentz- violating couplings to be 
unambiguously defined, they should be real traceless tensors with certain symmetry properties, such as H^^ — — 

It should be stressed that the theory remains ohserver-Lorentz-invariant, since the coupling constants are tensors 
under the boosts and rotations of the reference frame. This means that the Lagrangian is a scalar, and physics in 
different reference frames is equivalent though not the same, because the values of the components of the coupling 
constants depend on the choice of the frame. Nevertheless, these components in different frames also differ by the 
Lorentz transformation, hence, if one takes the latter one into account, physics in different reference frames becomes 
fully Lorentz-covariant. 
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In a fixed reference frame, however, the existence of the Lorentz violation affects, e.g., how the energy value of a 
free particle depends on its momentum, in such a way that its energy-momentum vector is no more a 4- vector under 
rotations and boosts of the momentum [active Lorentz transformations). In other words, extended electrodynamics, 
as well as the SME itself, is a theory with broken particle Lorentz invariance and preserved observer Lorentz invariance 

i- 

Most of the Lorentz- violating couplings are tightly constrained (see, e.g., Q and references therein); some of them 
can be excluded from the theory using certain unitary transformations, e.g., vanishes if we make a transformation 

-0 e-'^^^'V [I- 

In the present paper, we will study the contribution of the axial vector coupling yy'' into the electromagnetic (i.e., 
boson) Casimir effect between the two parallel conducting plates, i.e., within the theory with the Lagrangian 

C - -^F^.i^^'^ + VX^T^'^M-H^ + 'Ccs, (2.4) 
Ccs = ^7?^ 6^-/3^'^"" • (2-5) 

The CPT-odd Lorentz-violating correction Ccs entering the pure-photon sector of the theory has the form of the 
so-called Chern-Simons (CS) term [35|. Within the context of electrodynamics, it is also called after Carroll, Field, 
and Jackiw who have shown that this term can be induced by axions [3y| . It can also be a manifestation of a nonzero 
background torsion j37| . 

The theory with Lagrangian (12. 4p can be treated as a four-dimensional analogue of the so-called Maxwell-Chern- 
Simons (MCS) electrodynamics [35], which is usually considered in three dimensions and has the following Lagrangian: 

/:(2+l)D MCS QED - -^^M-i^''' + ^(^T'^^m " HV" + |e^a/3 (2.6) 

where is the three-dimensional Levi-Civita symbol, and /x, v^a^ (3 = 0, 1, 2. Due to this analogy, we will refer to 
the former theory as to the four- dimensional Maxwell- Chern-Simons electrodynamics. In (2-1-1) dimensions, the CS 
term does not violate Lorentz invariance, since ?7 is a pseudoscalar, and, when it is present, the free electromagnetic 
field F^i, satisfies the Klein-Gordon equation with a nonzero mass equal to 2|?7| f35|. However, the gauge invariance, 
which is typically absent for massive gauge fields, survives for nonzero 77, since the Chern-Simons term is changes by 
a total derivative under gauge transformations. 

The Chern-Simons term within (3-|-l)-dimensional electrodynamics has slightly different properties. It does also 
change by a total derivative under U{\) gauge transformations, and the action remains gauge-invariant (if appropriate 
boundary conditions are fulfilled): 

A,, ^ A^-d^,a{x), (2.7) 
^ ^ e'^^^f^V, (2.8) 

^ ^ f.-^eo.{x)^^ (2.9) 

C-d, (!|a(a;)eM-/3j^^^^ . (2.io) 
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However, first, the coupling constant 77^ is a pseudovector, and thus violates Lorentz invariance. In particular, it can 
cause a spatial anisotropy of the theory. Second, if r/^ 7^ 0, the theory can become unstable [J, H^. We will discuss 
this issue in Sec. IIVI and now we will consider the 4-dimensional MCS electrodynamics in more detail. 
Varying the action A — J Cd'^x with respect to A^i, i/j, and f/i, we obtain the equations of motion 

277^^^'' = f, (2.11) 
(7^(*9^-eA^)-™)V' - 0, (2.12) 
i^{^''ii1)^ + eA^)+m) = 0, (2.13) 

where j"^ = eipj'^ip is the fermion current, and the left arrow over 9^ indicates that it acts upon -0, i-e., to the left. 

As it can be seen from (12. 4p . the potential-current interaction has the conventional form —j'^A^ typical for the 
Maxwell electrodynamics. We therefore can resort to the pure-photon sector of the MCS electrodynamics for our 
consideration of the boson Casimir effect. The corresponding canonical energy-momentum tensor, which follows from 
the pure-gauge terms in the Lagrangian (|2.4p . reads 

T^T = d,A\F^^ + rfe^p^xA^') - v^^C, (2.14) 
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where r]fj.i, = diag(l, — 1, — 1, — 1) is the Minkowsky spacetime metric. Unlike the conventional (77^ = 0) case, the 
canonical energy- momentum tensor cannot be ultimately symmetrized so we present it in the following form: 

= 0,,^d^E[^,],, (2.15) 
S[am]„ = A.iF^x-eo^Pt.xV^A'^), (2-16) 
= F,xF\ + \rj^.F^pF^^ + re^p^pA^F^P [s^r^.x + l^H^li^^. - 5%.)) , (2.17) 

where Sj^^j^ — — S[^(;^]^. 6^^ becomes symmetric and gauge-invariant in the rj^ — Q case. 

Like in the conventional electrodynamics, the gauge freedom (|2.7p allows us to choose the transversal gauge (see 
Ref. [H) 

d^Ai" = 0. (2.18) 

However, in the rj^ ^ case, it becomes incompatible with the axial gauge n^A^^ = (n^ is an arbitrary constant 
4- vector) Q . The only axial gauge that can be fixed together with (|2.18p is 

r?^^'' = 0. (2.19) 

A specific feature of the MCS electrodynamics is that both T^|^" and even 6'^^ are gauge-dependent [l|. However, 
the Casimir force, which is determined through the normal-normal component of T^jlJ;" integrated over the plates, 
turns out to be gauge-invariant, for the special choice of rj^ we are interested in. We will prove this statement in the 
following subsection. 



B. Between the conducting plates 

Here and further on in the paper, we will assume that 

r;'^ = (r,,0). (2.20) 

We will consider the MCS electrodynamics of the free gauge field (without sources) between the two parallel plates 
with infinite conductivity. Let us choose the coordinate system in such a way that the plates are orthogonal to the 
z direction (i.e., to the 3rd axis, in accordance with — (t,x,y,z)) and correspond to z = ±a, where D = 2a 
is the distance between the plates. We will also assume that the plates have the form of large squares denoted 
P± = {—L/2 < x,y < L/2,z = ±a}, with linear dimensions L — > 00, and that the potential Afj,{x) satisfies the 
periodic boundary conditions, 

Af_,{t, x + L,y,z) = A^{t, x,y + L,z)^ A^{t, x, y, z). (2.21) 

It can be easily shown, using the method analogous to that adopted in the Maxwell electrodynamics 38], that the 
perfect plate conductivity implies the same boundary conditions for the electromagnetic field, as they are in the 
Maxwell electrodynamics, namely, 

E, = Ey = 0, xeP±, (2.22) 

where = Foi, i = 1,2,3, is the electric field strength (we will use the notation E = {E^, Ey, E^) = (E^ , E^ , E^), 
Ei = —E^ for three-dimensional vectors and denote their indices with latin letters i,j, k, ...). The Bianchi identity 

d^F"-' = (2.23) 

implies that 

doH = -TotE, (2.24) 

divH = 0, (2.25) 

where H is the magnetic field strength, iJ' = —■^eijkFjk ■ First of these equations, taken on the plates, gives 

doH, = dyE, - d,Ey = 0, xeP±. (2.26) 
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For the calculation of the Casimir effect, we can not consider the states in which the magnetic field has the static z- 
component that does not vanish on the plates (the corresponding one-photon eigenstates should have zero frequencies, 
but they, as we will show in (|2.45p . do not contribute to the Casimir energy). Another thing illustrating this is the fact 
that zero-frequency Fourier components of E and K are decoupled in the free MCS equations, so the zero-frequency 
component of the magnetic field is not constrained with the boundary conditions (|2.22[) . This means that this mode 
does not "feel" the plates and therefore does not contribute to the Casimir effect. As a result, we can assume that 

i7;,=0, a;eP±. (2.27) 

Now it is seen that, for our choice of rf ^ the change of the pure-gauge action between the plates under the gauge 
transformation (|2.7p vanishes: 

+00 



SA = J d^xSC = —^ J dt J dxdydz di (aeijkFjk) 



-00 V 

+OC +L/2 

dt J dxdya{t,x,y,z)H,{t,x,y,z)\lzt''a^O, (2.28) 

00 -L/2 



where V — {\x\, \y\ < \z\ < a}, i.e., the boundary conditions we chose do not violate gauge invariance. Of course, 
in the above expression, the function a{x) should be taken y-periodic as the potential Afj_{x) itself. 

The expression for the energy- momentum tensor becomes simpler when rj'^ contains only its timelike component. 
Of interest are the two following components: 

To'o^" = ^{E^ + H^)-'nA-H + E'diA°, (2.29) 

TsT = \{El + El-El + Hl + Hl-Hl) + ^Ti^^ (2.30) 

ATg^a^" = d^{A^F^'^)+7^{e:i,jA'dzAi + A- H ^2A^H'^). (2.31) 

The force acting upon the plates is a result of integration of T^f^ over them, i.e., over x,y E [—L/2, L/2], and 
subsequent averaging over the electromagnetic vacuum. To show that this force is gauge- invariant, let us prove that 
AT3I" vanishes after this integration and averaging. Indeed, one can easily show that 

J dxdyAT^I"" ^do J dxdyA^F°^ + J dxdy I ^ d^{A^F'^) ~ r]{di{A'^A^) - d2{A^A^)) j . (2.32) 

The vacuum expectation value of the first term, which is a time derivative, is zero, while the second one vanishes after 
the integration over x, y due to the periodicity conditions (|2.2ip . 

Hereby, we have proved that the Casimir force is gauge-invariant. To find the key expression for it, we can now 
resort to the axial-transversal gauge (|2.19p , (|2.18p , which, for our choice of 77^ , takes the following form: 

A° = 0, (2.33) 
divA ^ 0. (2.34) 

Within this gauge, the zero-zero component of the energy-momentum tensor reads 

= ^00 = liE^+ H^) H. (2.35) 

The equations of motion, in turn, take the form 

aA^2r]TotA. (2.36) 

Within this gauge, consider a complete set of one-particle electromagnetic field eigenstates. These states can be 
described by potential functions A = A„(a;)e^*'^"* satisfying the transversality condition (|2.34p and the equations of 
motion (I2.36p . which take the following form for these eigcnfunctions: 

{iol + V2)A„ = -277 rot A„, x e V. (2.37) 
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These functions should also be periodic in the y directions (|2.2ip and satisfy the boundary conditions (|2.22p . which 
imply that 

(A„), = (A„)j, = 0, a;eP±, (2.38) 

since we are not interested in the solutions with a;„ = which do not contribute to the Casimir effect. Finally, the 
eigenstates must be normalized in such a way that 

/ d'xAl{x)-A,Ax)^5^,n', (2.39) 
Jv 

where n, n' are the full sets of quantum numbers (excluding the energy sign quantum number) , which run through a 
discrete set since we assume the linear size of the plates L to be finite. 

As it can be seen from (12.371) . the equations of motion are time-reversible, since for every solution of the form 
A„(a;)e~"^"*, they also have a solution A„(a;)e''"''^"*. This comes from the fact that -q^ is a pseudovector, so its 
timelike component ?/ is invariant under time reversal. 

At the tree level, the electromagnetic potential between the plates can be quantized in such a way that 

A{x, t) = y] {A„(a;)e-'^-"*a„ + A*Jx)e^'^-*ai}. (2.40) 

Operators ajj/a„ commute in the usual way, 

[a„,a|j,] = 6n,n', [an,an'] = [ai,al,] = 0, (2-41) 

and create/destroy one photon in the mode with quantiun number n. In Sec. IIV[ we will show that, for sufficiently 
small |77|a, the imaginary-energy solutions (tachyon modes) form a zero-measure set, so the processes involving them 
are very rare, and the quantization (j2.40p is correct. 

To find the Casimir force /casimir acting upon a unit plate square (which is gauge- invariant, as it was shown above), 
we will differentiate the vacuum energy -Evac with respect to the distance between the plates. 



/casimir j- 2 7-» ' (2.42) 



1 dE,, 
dD 

Svac = (0| j d'xT^^^{x)\G) . (2.43) 

V 

Here, as usual, the negative sign of the force corresponds to the attraction of the plates. Let us show that the vacuum 
energy is half the sum over all positive one-particle eigenstate energies (as it is in the Maxwell electrodynamics). Note 
that for every solution A{x) of the equations of motion (|2.36p . due to the transversality (|2.34p of the potential, 

1 _ ^^^^ ^ A(Vx(VxA)) + V(Ax(VxA)) _ 

AjViVA) - V^A) + VjA X (V X A)) -AggA + V(A x (V x A)) 
= 2 i]AH = ^ , (2.44) 

and, after the integration over V , the last term, which is a spatial divergence, vanishes due to the boundary conditions 
(1135]). This implies that 

Sxn^ix) 10) = (01 j d^x (^^l±Ml r,AH^ |o> = (0| J d^x ^ |0) 

V V 

= E/ d?x^\{{~-iu:M{iuJnA*^)-A„{iio^fA*^)^Y. I ^'2; ^^^A„A: = i ^ a;„, (2.45) 

n y n y n 

where we have used the photon decomposition (|2.40p and the normahzation condition (|2.39p . Finally, 

/Casimir- - ^2 QJJ 2^ 2 ' ^ ' 
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The series under the derivative operator is obviously divergent, so it needs to be regularized and renormalized. One 
possible way to do this is to introduce the zeta function corresponding to the energy-squared operator that can 
be extracted from (|2.37p and has the eigenvalues equal to cj^: 

Kll = ^Ch^(-1/2), (2.47) 
C^.(.) = (2.48) 

n 

= -V^ - 27yrot. (2.49) 

The value C^2(— 1/2) is the analytical continuation of the series (|2.48p . which is convergent for large Res, s E C. 

Another possible renormalization method involves subtraction of an expression of the form (C^^) + C'-^^D), where 
are infinite constants, from the vacuum energy. C^^-* is obviously nonphysical, while C'^^^ equals the force acting 
upon the external side of the plate, i.e., the Casimir force in a semispace, so subtracting it is also physically motivated 
(see page [17]). Zeta function regularization automatically renormalizes the vacuum energy. 

We will use these two techniques in sections |Vl IVIl respectively, to find the corrections to the electromagnetic 
Casimir force due to the existence of a nonzero rj coupling. 



III. ELECTROMAGNETIC FIELD MODES BETWEEN THE PLATES 



In this section, we obtain the expressions for the solutions A„ {x) of vacuum Maxwell-Chern-Simons secular equation 
(I2.37P with Ty'' = (77, 0), i.e., solve the free photon eigenstate problem. 

The translational invariance along the x, y directions allows us to search for the solutions having the form 

An{x)=Ne'^''f{z), k = {kx,ky,0), (3.1) 

where N is the normalization constant. Due to the invariance with respect to rotations on the z axis, we can choose 
the cartesian coordinate system so that ky = 0, kx — k > 0. Of course, the periodicity conditions (|2.2ip with L 00, 
which are not invariant under these rotations, do not affect the form of the solutions, but only the normalization 
coefficient N and the possible values of k, so that k^.y = ^nx,y, nx,y € Z. 
The gauge-fixing condition ()2.34p takes the form 

kfx = idjz. (3.2) 

On the boundary (i.e., plates, z = ±a), Ax^y = 0, hence 

/,,,(±a) = 0, (3.3) 
djz{±a) = 0. (3.4) 

The projections of the MCS equations of motion (|2.37p onto the x, y, z axes read, respectively, 

(fc2 _ _ a^)/, ^ -2fjdzfy, (3.5) 

{k'-ujl-d'z)fy = 2ri{-ikfz + dzfx), (3.6) 

{e-ul-dl)fz = 2r,-ikfy, (3.7) 

where a;„ is the energy corresponding to the eigenstate. Note that differentiation of (|3.7p with respect to z and 
subsequent substitution dzfz = —ikfx (according to (|3.2p ) gives (|3.5p . so we can discard the latter equation from our 
consideration. 

Using p.2p . substitute f for dzfx in (|3.6p . then we come to the system of equations on fy, fz- 

{e-u:l-dl)kfy - -2zri{k^-dl)fz, (3.8) 
{e-ujl-dl)fz = 2ir,kfy. (3.9) 

These equations, together with the boundary and gauge-fixing conditions, possess a special kind of discrete symmetry 
(we will call it the z-parity and denote 11): 

(2kk\ 
kj j Aj{x,y,-z), k = {-idx:-idy,0), (3.10) 

Ufy,z{^) ^ fvA-^), Ufx{z)^-fx{-z). (3.11) 
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The transformation law of the x-component, compared to that of the j/, z-components, is in accordance with (j3.2p . 
Let A{x) (and f{z), in turn) be the eigenfunction of the z-parity with the eigenvalue 11, which is obviously +1 or 
— 1, since iP — 1. 

The solution f{z) corresponding to definite 11 and fc is a superposition of the two solutions f^^'^\z) in the boundless 
space (i.e., without the boundary conditions) with the two polarizations A = ±1. As it is seen from p.Sp and (|3.9p . 
we can search for the y, ^-components of these solutions in the form 

/(n/)(z) cx ^ni^z), (3.12) 

= |sin>.z, n = -l, (3-13) 
where ><• is a complex constant. Then the equations (jS.Sp . (|3.9p are easily solved, and we obtain 

-^AcosJ^^(-^^)' ^'"^±1, (3.14) 




= -ryA+Vc^T^, (3.15) 
Kl~-k^, (3.16) 

cosOx = sin^A-^. (3.17) 

As seen, the sign of x\ does not affect the form of the solution, except for its sign. Note that the solutions with 
complex >fA, which are unbounded for \z\ oo, should not be discarded from consideration when finding the solutions 
between the plates, where z is finite. In the next section (see pagellOp. we will show that >cx can be either real or pure 
imaginary for these solutions. 

Now consider a solution of (|3.8[) . (|3.9p with definite 11 = ±1 satisfying the boundary conditions (|3.3p . (|3.4p 



fyAz)^ E CA/("/n^)- (3.18) 

A=±l 

Since the boundary conditions at z = +a and z = —a are equivalent for the functions with definite z-parity 11, it is 
enough to consider these conditions only at the point z = +a. In the following expressions, we denote ^± = for 
the symbols depending on the polarization A, such as C\, k\, 6\, etc. Then the boundary conditions for fy and fz 
take the form: 

C+iy9n(>^+a) + C'_(pn(>?-a) =0, (3.19) 

C+ • (— i) cos0+(— ^>^+ly9_^(>^+a)) + C*- • * cos(?_(— n><_(p_n(>^-a)) =0, (3.20) 

where we have used the identity 

dzifn{x^z) — —Tlxip^Yiixz). (3.21) 



The boundary conditions for are satisfied in turn, due to (|3.2p . The existence of nontrivial solutions for C± requires 
that the determinant of the system of equations ()3.19p . (I3.20p vanishes, then 

gui^Jn) = ipn{}<+a)(p-nix-a,) sinO^ + (^n(>*-a)</'-n(^+a) sin6'+ = 0. (3.22) 

This equation implicitly determines the spectrum for every fc > and 11 = ±1. When it is solved, the solutions for 
C± read (up to a nonzero multiplicative constant) 

C± = ±(pn{>c^a). (3.23) 

Now, to determine the final expressions for /(z), we use the solutions without the plates p.l4p and then Eq. p.2p to 
find fx{z). As a result, we obtain 

/(z) — —Ex ■ n(sin6'+(pn(>f-a)<;5-n(>f+2) + sin6'_(/9n(>f+a)<y5-n(>f-2^)) + 

+ey ■ {ipn{>i-a)ipn{>i+z) - ipn{^+a)ipn{x-z))~ (3.24) 
—Bz ■ i{cos9+ipu{x-a)ipui>c^z) + cos6^ipYiix^a)^nix-z)), 
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where, for arbitrary fc = (k^, ky, 0), 

k 1 

and Bz is a normal to the plates pointing in the direction of positive z, i.e., = Vz. The full set of quantum numbers 
is n = {kx, ky,Il,m), where k^ = 2^^, ky = n^.y G Z, 11 = ±1, and the last quantum number m = 0,1,2,3,... 

enumerates the solutions of the spectrum equation p.22p . 

When 7] = 0, i.e., in the Maxwell case, (|3.22p gives the well-known spectrum 



K)m = (^)', m = 0,1,2,3. (3.27) 



The solution for the potential, corresponding to m = and 11 = —1 is trivial in this case, the others are physical. 
However, the energy eigenvalue with m = 0,11 = +1 does not depend on a and hence does not contribute to the 
Casimir energy. 

IV. GENERAL NOTES ABOUT THE ENERGY SPECTRUM 

In this section, we will show that the Maxwell-Chern-Simons photon energy spectrum between the plates is such 
that the energy-squared eigenvalues are real, > —rf', and, when |r/| < the imaginary-energy solutions have 
cij,^j = —rf and form a zero-measure set in the quantum number space, i.e., almost no tachyon modes are present. 

Consider a space J- of potential functions A(a;), defined in a spatial domain V between the plates, which are 
transversal l|2.34p and a;, y-periodic (|2.2ip . The scalar product and the norm in T are defined in the usual way. 



{A'\A)^ d^x{A{x)rA{x)^ IIAII ^ yp^). (4.1) 



It can be easily seen that the energy-squared operator H"^ (see Eq. (|2.49p ) acts inside of since VA G !F : 
Ayv{H'^ A{x)) = . Denote the orthonormalized eigenstates of H'^ in T (which were discussed in the previous 
section) as |An), so that \ An) = \ An)- The corresponding eigenfunctions An[x) should also satisfy the bound- 
ary conditions (|2.38p on the plates P±, which, in turn, imply that dz{An)z — divA„ — dx{An)x — dy{An)y = on 
them. Let us prove the following assertion. 

Assertion 1. The eigenvalues uj'^ are real and uj^ > —rf'. 

Proof. Let us denote = (^^| rot |A„) and = (AJ (-V^) |A„), then = (A„| H'^ |v4„) ^ (Ai|-V2-277rot \An) = 
— 2ri^s- First, integrating by parts, we obtain 

= j{A^)*{~djd,{A^),)d^x^ j \d,{A^\\''d^x~ j{A^)*dMn\dSz. P = P+UP^, (4.2) 

V VP 

where the area element dSz — idS — ±\dxdy\ on P±. Indeed, due to the periodicity, the integration by parts over x, y 
does not give the surface term analogous to the last one in the above expression. The latter term, however, vanishes 
due to the boundary conditions on P±, namely, {An)x = {^7i)y — 0, dz{An)z — 0. Thus, 

^p = J Id.iAnhfd^'x e [0,-Hoo), (4.3) 

V 

i.e., real and nonnegative. Making the analogous integration by parts in the expression for ^s, we obtain 

6 = y" e^Jk{An)*^J{A„)kd^x = J {-dj{A,,)*)eijk{An)kd^x + J ei3kiAn)*{A^)kdSz, (4.4) 

V V P 
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and e,3fc(A„)*(A„)fc = (A„);(A„), - (A„)*(A„)j, = on P±. Finally, 

6 - / {~dj{A„)*)e,,kiAn)kd^x = [ (A„)feefe,.(5,(A„),)* = 4% (4-5) 



meaning that is real. Thus we have proved that ui^ = — 2?7^s is real. 

Now let us show that > — 7y^. Consider the two states \A') = \An) and \A) = rot|A„). Note that \A) £ T. 
Then, according to the Cauchy-Bunyakovsky-Schwarz inequality [33|, | (A'|A) P < i.e., 

il = I (A„| rot |A„) |2 < II rot |A„) ||2 = ^ (A„)*e„,, 5 ,e,,fcab(A„)fcd3x = ^ (AO* (9 ,5,^9, - ^9 fc9,)(A„)fcd33,. (4.6) 

V V 

Now, the integration of the second term by parts makes d operators act to the right, and we obtain 

J \da{A^)k\^d^x + J{A^)*dkd,{A^)kd^x- JiAr,)*d,{A^),dS,, (4.7) 

V V P 

where the second term vanishes due to the transversality, while the third due to the boundary conditions. Finally, 
comparing with (j4.3|) . we can conclude that 

e < ^p- (4.8) 

Hence, the energy-squared eigenvalue 

^I^^P- 2^76 >^P- 2\v\VCp ^iV^- hi)' 'V^> ~V^- (4-9) 

□ 

As seen, the assertion proved above is also valid in the infinite space (a 00). Then, since ojf^ + > 0, the two 
branches of the spatial momentum K\, which are defined in (|3.15p . are real, while x\ = ^ — can be either real 
or pure imaginary. Moreover, one can easily see that 

iol = K+K^, (4.10) 

i.e., the imaginary-energy solutions correspond to K± having opposite signs, and the real-energy solutions correspond 

to K± > 0. We can also note that cos 0a = k/K\ is real however it may not lie within the segment [—1,1], and 
sin6'A = x\/K\ is either real or imaginary (the latter case when | cos 0^1 > !)■ 

Further we show that, for sufficiently small finite a, the imaginary-energy states are negligibly few. 

2 

Assertion 2. For the solutions with negative z-parity 11 = —1, > . 

Proof. As it was shown in the previous section, we can choose the coordinate system in such a way that 

A„(a;) =iV/(z)e^'=% A: > 0. (4.11) 

Then, using (|4.3p . we obtain 

||A„||2 = \N\'L^ff, (4.12) 

^p ^ (A„| (-V^) K) = J\dMn)A^d'x=\N\^L^k^fr + \\djr), (4.13) 

V 

a a 

Wff ^ J\f{z)\^dz, \\djr ^ J \dj{z)\^dz. (4.14) 

— a —a 

Due to the boundary conditions (|3.3p . (13. 4|) . the negative 2;-parity fif{z) = —f{z), and the transversality p.2p . we 
can write f{z) as a Fourier series on the segment z G [—a, a]: 

„ , , i 7r(m — 1/2) 7r(m— l/2)z 
W - tY.- —a.^mcos^ 4.15 



m— 1 

00 



iy\^) ^ 2^ay,mSm , (4.16) 

m— 1 

/,(z) = ^a,,„sin^ (4.17) 
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Using the fact that 



irmz . 7rm z 



sin sin dz — / cos 



irmz irm z 
- cos 



dz = ahn 



. 7r(TO- 1/2)2 . 7r(TO'-l/2)z, 
sm sm dz — ao„ 



7r(m-l/2)z 7r(m'-l/2)z, , , , „o 

cos COS dz = ao„i m' , m,m =1,2,3,, 



(4.18) 
(4.19) 
(4.20) 



we find the expressions for the norms in terms of the a-coefficients: 

'7r(m-l/2) 



II/IP = aEll" 

m— 1 



y,m I 



1 



7r(TO- 1/2) 



(4.21) 
(4.22) 



Now note that \ ^ay^m\ > 



7r(m-l/2) 



> 



^|ctz,mP, then the foUowing inequality takes place: 



\djr>^\\fr, 



(4.23) 



which, due to (|4.12p . (|4.13p . leads to the lower constraint for in a normalized state |A„) with negative z-parity: 



(A„|(-V^)|A„) 

IU„||2 



\dj\ 
\\f¥ 



(4.24) 
□ 



Together with (|4.9p . this obviously leads to the statement. 
Assertion 3. When a < the negative z-parity eigenstates correspond to the eigenvalues > 0. 



The method we used above to place the constraint on the energy spectrum of the negative z-parity states does 
not apply to the states with positive z-parity. Indeed, the state with /(z) = (0,0, /z), fz = const, which satisfies 
the gauge and the boundary conditions, corresponds to = fc^, which can be arbitrary small nonnegative number. 
This is inconsistent with the method we have used in the proofs above. However, it is possible to derive a different 
spectrum constraint directly analyzing the properties of the equation (7_|_(a;^) = (see Eq. p.22p ). which implicitly 
determines the energy spectrum of positive z-parity states. 

Assertion 4. When a < 2j^, the only nontrivial solution with = —rj^ has the form f{z) = (0, 0, const) for k — \ri\. 

Proof. Note that when one replaces 77 with — r/, then K+ ^ and all quantities that depend on the polarization 
A — ±1 flip in the same way. The functions gYi[LO^) are symmetric with respect to this polarization flip, so they are 
also invariant under the change of the sign of 77, and the properties of the energy spectrum, in turn, do not depend 
on it. In view of this, we will assume that 77 > in our further speculations. 

For — —rf^ (minimum possible value), K± = ^r], >f+ = x- = yTy^ — fc^, and g+ = 0. However, we should also 
analyze the existence of nontrivial solutions for /(z), and not only for C± (see Eq. p.23p ). Then let us write the 
expressions for fy,z{z) for this case (defined in the previous section): 

fvA^) ^(^y.z cos (^^ if -k'^ z) , \ay\^ + \a,\^ >0. (4.25) 

For k > T], this solution is inconsistent with the boundary conditions (|3.3|) . (|3.4|) . When k = rj, the only solution is 
the one announced in the assertion. 



fy{z) = 0, fz{z) = const. 



(4.26) 
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This solution exists for any a > 0, so it makes a constant contribution to the vacuum energy, and thus does not 
affect the Casimir force. Moreover, the solutions of this type form a zero-measure set compared to the complete set of 
eigenf unctions, since they exist only for k — rj, while other solutions exist within a finite-measure set in the fe-space. 
This means that the contribution of the quantum processes involving these states is negligible. 
In the k < rj case, the boundary conditions are satisfied with the following two types of solutions: 

/(iHz)=0, fP(z) = const- cos k = Jrf - (nm/a)^ > 0, (4.27) 

" a 

(z) ^ const ■ cos ""^"^ ^ ^^^^"^ , fP(z)^0, k = Jrf - {7r{m - l/2)/a)^ > 0, (4.28) 

where m — 1, 2, 3, ... When a < the solutions are inconsistent with k > for any m — 1,2, 3,..., so there are no 
solutions of this type. Therefore, in this case, the only solution with ujf^ = —i]^ has the form (|4.26p . □ 

Assertion 5. When a < the solutions with —rf' < cj^j < do not exist. 
Proof. Let us introduce the two functions: 

5n(w^) = , , r = tan" x+a sin0+ + tan" x^a sin6'„. (4.29) 

ipu{>c+a)Lpii{x^a) 

Let us treat gn and gn as functions of Kj^ G [— ?7, -|-oo), with K_ = -I- 2?], ?7 > (both gn and gn do not depend 
on the signs of ?/ and oj). This interval is in the one-to-one correspondence with o;^ = K+K- € [—ff, -l-oo). Consider 
the domain with negative uS^ £ {~ri'^,0), which corresponds to e (— ry, 0). Here, >c± e (+ioo,0) U [0, 2?7) and, 
because we assume that rj < cos>c±z ^ for any z £ [—a, a], so the denominator in the expression (|4.29p for 
is nonzero. Thus, the equations (7+ = and 5+ = are equivalent in the domain X+ e (—77, 0), moreover, ^_|_ is a real 
continuous function of Kj^ E (—77,0) for any fc > 0. At the end points of this interval, 

5+|x+=-r,+ = 0, (4.30) 
5+k+^o_ ^ (4.31) 
Taking the derivative of with respect to wc obtain 

= V f 1 + tan^ + '-^^ cos^ 0,) . (4.32) 
a oK+ ^ \ x\a I 

^ \=±1 ^ A / 



Here, cos^ 9± = k'^/K^. > 0. On the other hand, if >cx is imaginary, then 

> 0. (4.34) 



l + tan^XAfl = 1 - tanh^ |>fAa| > 0, (4.33) 
ta,nx\a tanh \ x\a\ 



If >i\ is real, then < yi\a < 1r]a < 7r/2, hence tan>c\a > 0, and 

1 + tan^ >cxa > 1, (4.35) 

'-^^ > 0. (4.36) 
xxa 

Finally, we conclude that the whole expression ()4.32p is positive: 

dg+ 



> for K+ e (-77, 0) and fc > 0. (4.37) 



Taken together with (|4.30p . (|4.3ip . this proves that g^{u!'^) is positive within the interval e (— r/,0), and the 
corresponding one-photon eigenstates do not exist. □ 

Thus, we have come to the following conclusion. 
Assertion 6. When a < the only negative-uj"^ Maxwell- Chern- Simons one-photon energy states are the following: 

A^{x) = Ne^'^'^e,, (4.38) 
ivl = -77^ (4.39) 
fc = (fc,,fey,0), |fc| = hi, = (0,0,1). (4.40) 

These states form a zero-measure set in the k-space. All other states are non-tachyonic (i.e., correspond to ui^ — ^)- 
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V. THE CASIMIR ENERGY: ZETA FUNCTION REGULARIZATION 



In this section, we will find the corrections to the Casimir force and the Casimir vacuum energy calculating the 
leading 77-correction to the zeta function defined in (|2.48p . in the case |7y|a ^ 1. As it was shown in the previous 
section, in this case, the solutions with imaginary energies form a zero-measure set in the fc-space, then the reduced 
zeta function, defined per unit plate square, 



1 



1 



12 



kdk 



2Tr ^ ^ 

n=±i m 



('^fe.n.rn) 



(5.1) 



does not include the contribution from these states which vanishes divided by L"^ ^ 00. On the other hand, the 
contributions from all real-frequency states are unambiguous, since (w^)"" = e~^^°^^" when oj^ > 0. Here, n — 
{kx, ky, n, m) is a full set of quantum numbers (see Sec. IIIip . The series (|5.ip typically converges for sufhcicntly large 
Res and thus can be analytically continued to s G C, in particular, to the point s = —1/2, which corresponds to the 
renormalized sum over the frequencies ujk,u,m- 

As mentioned in the previous section (see page fTTj) . the frequencies Un do not depend on the sign of ry, then ({s) 
is also an even function of 77, or is, in other words, a function of i]^. Then the derivative of the zeta function with 
respect to i]^ reads 



d 



d{rf- 



-as) = 



kdk ^ — ^ ^ — ^ 



-2s 



n=±i m 



(5.2) 



Let us now find the derivatives doJkji,ml d{rj^) near the Maxwell solutions (i.e., those in the 77 = case, see Eq. p.27p ) 
for m = 1, 2, 3, ... explicitly showing the dependencies of the functions on rj in the following expression: 



^5n(t^„(?? ),?y)|r,=o 



\ drj 



d'^gn{ujl{rf),ri)duJn{ri^) dgYi{ujl{Tf),'q)d^uJn{rf) 



dujndrj 



drj 



dujn 



(5.3) 



j)=0 



Since ujn is an even function of 77, its first derivative dujn/drj vanishes at 77 = 0, thus the second and the third terms 
vanish, and we obtain 



7( = 



2 



drj"^ 



ri=Q 



2 dgii{i^l,ri)/dujr. 



(5.4) 



J7=0 



Now, to find the second 77-derivative, let us treat gn as a function of K± rather than a function of uj and r/. Then, 
according to (|3.15p . 



dgn 
du! 



ri=0 



V = 
d 


-K+ 


2 

dK+ d 


drj 




drj dK+ 






\{ ' - 


ri=0 




\dK+ 



^JK^, gniio^ rj) = Gn{K+,K-), 



dK^ d 



2K, 



d 



2K_ 



d 



drj dK^ K+ 
d V I f d 



dK- 



dK 
d 



dK^ ' 



K+ \dK+ dK- 



K+=K-=u 



d'^gu{(^,(^) 1 dgn(^,^) 
cLu'^ UJ dw 

dgn(^,^) 
duj 



-4 



d'Gu 



dK+dK. 



K,=K_ 



(5.5) 
(5.6) 



(5.7) 
(5.8) 
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Taken at the zeros ujm of gn in the ?/ = case (see Eq. p. 271) ). these derivatives read 



drf 
dgn 



(-l)"2a 



17=0 



i-(i-2(-irn)-^ , 



^^ + ^M' ^^^~2^' 1,2,3,. 



Substituted into ()5.4I1 . this results in the following expression for the frequency shift: 



d{rf) 



1 



))=0 



i)=0 



2lom 
1 



i-(i-2(-i)"n)-^ 



1 



(5.9) 
(5.10) 
(5.11) 

(5.12) 
(5.13) 



Remembering that in the Maxwell case, the frequencies com do not depend on the parity 11, we substitute this result 
into ()5.2|) and, after the summation over 11 = ±1, obtain: 



d 



dm 



_ r kdk ^ 2s / 4a2A;2 

^'\=o ~ J "2^£'i(P + (W2a)2)^+^ 

= y ^ ^ (P + ( W2a)^)- = 

„ m— 1 ^ \ / / / 



4a^ 

-(/(s + l,0)-—/(s + 1,2)1, (5.14) 



(5.15) 



00 

/oo 



(k^ + (7rm/2a)^) \m 



(5.16) 



These integrals can be calculated using the integral representation of the Euler gamma function, the Gauss integral 
identity, and the series representation of the Riemann zcta function Cr(s) [40| . 



(fc2 + (7rTO/2a)2)' 



r(s) 



r § + 1 

, Rea>-2, i>0. 



00 ^ 

V — 



Cr(2s-2), Res > 3/2. 



(5.17) 

(5.18) 
(5.19) 



Let us evaluate the integral I{s, a) using these three identities referred to above the equality signs where they are 
used: 



J(s, a) 



sTrt 1 



/ kdk 



Vis) 

^ ' m— 1 

isTTst r(f + 1) 
2r(s) 

r(f + i)r(s- f ~ 1) (2a 

2V{s) 



1 ^ 1 



r(s) ^ TO° 

^ ^ m=l 



Enj r(f + i)r(s - f - 1) ^ i /2a 



2r(s) 



E 





2S-Q-2 



2S-Q-2 



Cr(2s-2). 



(5.20) 



These transformations are valid (i.e., all integrals are convergent) when Rea > —2, Res > max(|, 1 + Re ^), i.e., for 
sufficiently large Re s for any fixed a, such that Re a > —2. 
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Then, for Res > |, both the zeta functfon <^{s) and its ry^-derivative at 77 = are convergent, 



2 ' 

1 /2a^ 

2tt{s - 1) V tt 

d 



C(s)|,=o = ( - ) Cr(2s - 2), (5.21) 



d{7f) 



s-2 /'2a^^" 



^=0 27r(s - 1) V I" 



— Cr(2s), (5.22) 



as) = ^^(^ (7) ' (cr(2s - 2) + (s - 2) (^^) ' Cr(2.) + o{ifyj , (5.23) 

where Z? = 2a is the distance between the plates, as it was mentioned before. 

One should also analyze the contribution of the "quasi-zero" modes into the Casimir effect. These modes formally 
correspond to m = and can be left aside in the ij — case (see comments after Eq. p. 271) ). It is easily seen 
from (j5.13p that this formula does not not have its asymptotical meaning in the m = case, since >cm = 0. Proper 
asymptotical solution of the equation (|3.22p near uj = k (i.e., for m — 0) gives the following expressions for the energies 
^fc.n,o of the "quasi-zero" modes with parities 11 = ±1: 

^fc,+i,o = fc fl - ^V^a^ + 0{rj^)) , k » r,, (5.24) 



2 

a;fc,_i,0-fc±7?-|^+O(r?=^), A: > 77. (5.25) 



The odd (11 = —1) solutions do not contribute to the Casimir force within the second order in 77 since they do not 
depend on a. The even (11 = +1) solutions do depend on it, and the corresponding ry-dependent correction to the 
zeta function C,{s)^ taken for k ^ r/ solutions, reads 



2 / kdk —2s / 2 2\ 2sTfo? f dk s rfa? 



which has the order 0(77^) when s —1/2. The corrections to the zeta function due to the domain k S [0, 0(77)] have 
the order o(j]^) due to the smallness of this domain. Thus the "quasi-zero" modes do not contribute to the Casimir 
force, within the second-order accuracy in rj. The strict demonstration of this fact will be presented directly in the 
next section, where the exact integral expression for the Casimir energy will be derived. Now we limit ourselves to a 
qualitative discussion of this issue. 

Using expression l|5.23p . we finally obtain the correction to the Casimir energy and force. 



-Evac _ C(-l/2) _ , o^ ^i^D) 



^2 



Cr(-3) - ^^Cr(-1) + oiivD) 



2 V 27r2 

oiivD f) ] , (5.27) 



This correction is relevant when the dimensionless parameter \rj\D <C 1, i.e., at relatively short distances. It should 
be stressed that the zeta function regularization automatically gave us the renormalized result for the Casimir energy, 
within the leading order in rj. 



VI. THE CASIMIR ENERGY: SERIES SUMMATION VIA THE RESIDUE THEOREM 

In this section, we will make a renormalization of the vacuum energy series (12.451) using the direct summation of the 
one-particle eigenstate energies implicitly defined as the zeros of the functions .gn(w^) (see Eq. (|3.15l) _). The explicit 
expression for the sum over the energy eigenvalues Wfe,n,m can be derived using the residue theorem [34| . 

First, as it was mentioned in Sec. llVi the properties of the spectrum do not depend on the sign of 77, so we assume 
that 7; > 0. Then, it was also shown that the energy eigenvalues are either real or pure imaginary, and all of them 
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correspond to G [—77,00). The real eigenvalues correspond to (E [0,cx)), so the real part of the vacuum energy 
can be presented in the following form: 



1 ^ 1 /" kdk „ v-^ „ 

n=±i 

1 

Sn{D) — ^k,m,n, 



(6.1) 
(6.2) 



where means a summation over those m, for which uJk,m,n G (0, 00), i.e., the corresponding e (0,00). Let us 
apply a regularization and redefine Su as 



(6.3) 



where A —^ +00 is the cutoff parameter. This sum is convergent when A is finite and can be represented as a complex 
plane contour integral. Recall the functions introduced in section HVl instead of f/n(w^): 



ffn(w^) 



5n(w2) 



= tan'^ x+a sin 9+ + tan^ x-a sin 0- . 



(6.4) 



These functions can be treated as functions of E C (for fixed fc > 0, as we will further assume). Al 



ReK. 




FIG. 1: To the transformation of the sum into a complex plane integral. The crosses indicate the zeros of the function gn, 
which correspond to the real energy eigenvalues. 



though >c± = ^ K\ — have two branches (namely, they are defined with an arbitrary sign) , the expressions 
tan'^ >!r±a sin 6*^ = ■^x±tan^ x±a have only one branch, moreover, functions have only pole singularities, i.e., 

they are meromorphic functions [s^ . Note that the original functions gn were multivalued complex functions of , 
so they were incompatible with the residue theorem. 

Since functions (pn{x±) are bounded in any finite circle {l-ft'+l < M}, M > 0, the set of zeros of gu lies within the 
set of zeros of gu- It implies, in particular, that the zeros of the functions also correspond to real values of K^, as 
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it was for the zeros of gu- Moreover, it can be explicitly shown that the difference between the two abovementioned 
sets contains only those points for which either G {0, ±fc}, or K_ G {0,zLk}, or (pYi{x+a) = (pYi{>C-a) = 0. The 
first two cases correspond to the energy eigenvalues that do not depend on a and therefore can be excluded from the 
Casimir energy series. On the other hand, the solutions of the equation (^n(>^+a) — (/3n(>f-a) — do not exist for 
almost all k (all but no more than a countable set), so the summation over these additional energy eigenvalues gives 
a vanishing contribution to the vacuum energy, after the integration over k. Thus it is convenient to think of Su as 
of the sum over all zeros of the function ^n(w^) corresponding to G (0, oo), regularized in the way shown in (j6.3p . 

Now we can apply the residue theorem to the domain A = {ReKj^. > 0, \K^\ < A} (see FigllJ. Denote the contour 
enclosing this domain as C = Cjm + C\, where Cim corresponds to the segment of the imaginary axis, and Ca is the 
semicircle with the radius A, as shown in Fig[TJ The direction of this contour is chosen "mathematically-positive" , 
i.e., counterclockwise. Then the residue theorem implies that 



dK+ 
27ri 



(n)^„-a;(/f<"')/VfcA 



Res 



(6.5) 



where Ld{K+) = ^JK^K- = ^ + 2if) is the analytical function in A. K 
lying within A, while k'^F'^ correspond to the poles of the numerator, i.e., dgu/dK-^ 



(n) 



correspond to the zeros of gn 
The first sum in the right side 



of (j6.5p can be replaced by SyyD when A — > oo, since we have shown in the previous paragraph that we can make a 
summation over the zeros of the function gu instead of gu- Let us substitute the functions and their derivatives 



1 d~gn 
adK+ 



- E n(i 



A=±l 



tan^'^ x\a) 



tan'^ >c\a 
XX a 



cos^ 9\ 



(6.6) 



into (|6?5]): 



1 / ^ ),-^(K,)/VkA d~9n/dK+ _ 1 y ^(^(rt)^-.(Ki 
la J 2Tn gji 2 ^ 

c KY' 



X Res 



E (n( 



1 + tan™ ^^a) + ^'^"^ cos^ 6»a 



sin 9+ tan^ x+a + sin 9- tan^ >i-a 



-,K+=K_ 



(p) 



(6.7) 



One can see that the summation is performed over those values of for which either tan'^ x^a — > oo or tan^ x^a 
oo, and the residue is equal to the residue of the function \E'n(if+) = H + *'^"in ° ) ' "^^i^h does not have 

any singular points other than K^^"^ . Then the sum over the poles K'f^ in the right side of (|6.7p can be also presented 
as an integral over the contour C, with an integrand containing '^yi{K^). Then the two integrals, taken together, give 
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2a 



'2TTi 
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n f dK 
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,-^(K+)/VfeA 



d~gn/dK, 



a*n(ivr4 
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A=±l 



27ri 



u;{K+)e- 



.(K+)/VlA£:n 



1 — tan'^ xxa ta,n^ x^\a 



sin 9\ 



n tan'^ K\acos^ 9^ 



sm w_A 



(6.8) 
(6.9) 



Now we need to renormalize the integrals we have obtained, i.e., subtract their A — > oo divergencies, which can be 
extracted analyzing the case a — > oo, i.e., infinitely distant plates. Two types of divergent counterterms should be 
subtracted, namely, in terms of ^n. 



(6.10) 
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(1 2) 

where Cfj ' are the functions of k, r], and the regulator A, but they should not depend on the distance D between 
the plates. The term with gives a constant contribution to the vacuum energy, so it is non-physical. The other 

(2) 

term, which contains , gives a constant contribution to the "average vacuum energy density", that is, the vacuum 
energy divided by L^D. It is obvious that, for large D, this "vacuum energy density" is equal to the force acting on 
a plate in a semispace. Thus subtracting it means the subtraction of the force acting upon the open side of a plate, 
which is physically reasonable. 

First, let us analyze the asymptotic behavior of the integral over the semicircle Ca for A ^ cx). When = Ae*'', 
\d\ < tt/2 — e, e > 0, the integrand is exponentially small due to the factor g-^/^^i^ {uj{K+) is a one-valued function 
in A, with Rew > 0). On the other hand, when is near ±7r/2 so that Reuj ^ VkA, we can use another set of 
approximations: 

u;,K±,K± = 0{A), (6.11) 
sin6l± = 1 + 0(1^2), (6.12) 
cos6i± = 0{l/A), (6.13) 
tann>f±a = sgni? + ©(e^^Aa-^^ ^q -^^^ 

where = + 2rj. Then the integrand 

u;(i^+)e-'^(^+)/^=ii =cj(i^+)(-4msgnt9 + 0(l/A2))e-"(^+)/^, (6.15) 

moreover, the terms entering this expression that depend on a (in other words, change of this expression if it is 
treated as a function of a) have the order ©(l/A^). Since (jj{Kj^) = 0{h) and the integration effectively includes the 
pieces of C\ of length 0{\fkK), where the exponential factor is not infinitely small, the integral over Ca has the order 
0{VkI^) while its parts that depend on a have the order 0{^Jk/ K) and vanish in the A ^ oo limit. The latter fact 

(2) 

holds true even for a ^ oo. Thus the integral over Ca has the form of the countertcrm containing Cj^ , except for 
the terms that vanish in the A — !■ oo limit. Hence, this integral is fully cancelled when renormalized. 

Second, let us consider the a ^ oo asymptotic of the integral over Cim to extract the counterterms. In this limit, 

tan'^ >f± a — > illcr, (T = sgn ImiC+, (6.16) 

where the difference between the left and the right sides is exponentially small in a, i.e., smaller than any negative 
power of a. Then the integrand takes the following asymptotical form: 

Hn f I 1 \ n fcos^e+ cos^e\ 

-7— = -iUa - — h - — — H — — : — — — \ + exp-small terms m a. (6.17) 

gu \smt/+ smt^_/ a[smO-^- + suit)-) \ >t^ J 

Here, we do not need to strictly specify what the exponential smallness in a means (e.g., write expressions of the 
form ©(e"^"' ^"^±1)), since our asymptotical treatment is only aimed at extracting the relevant integral representation 
of the counterterm, which should be compatible with (|6.10p . Namely, we see that the first two terms in the above 
expression contribute to C^^"* and respectively. Hence, we subtract these two terms, multiplied by Loe"^^^^^, 

from the integrand tJe~"/^'^*^Sn/5n fof finite a. Since the integrand over Ca vanishes after renormalization and 
setting A —^ oo, then, after this subtraction, we obtain the integral representation of the renormalized sum over the 
frequencies: 

z J Zm sm 0+ tan >ir+a + sm 0- tan K^a 



^ren 



o 4- n j, n sin6'A / sin6iA \ HAA' cos^ 6'a' sin6l_A' 1 

2 - > tan xxa i tan K-xa^- zHcr 1 + -— — - > — — - — } . (6.19) 

I sm0_A V smO^xJ ^^^^^ (sm6'+ -|- smt'_)>fA'a J 



A=±l 



One can easily show that the integrand is exponentially small for | Im_ftr^_| 3> »7, 1/a, even without the regularizing 
factor e^'^/^/*^. Moreover, the integral is convergent near = 0. Now we are able to set the regulator A to 
infinity, then the smooth cutoff factor ^-^^I^J^ disappears. The integration is performed over the imaginary axis 
Kj^ G Ciin = (-fioo, — zcxd), so it is natural to make a reparametrization making Kj^ real (namely, Kj^ — > 
dK^/2'Ki —dK+/2'K, G R). Moreover, looking at the expression (|6.ip for the vacuum energy and the expression 
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for S'n", we see that we can also make a reparametrization of the integration variables making them dimensionless, 
such as ka and K^a. Making these two types of reparametrizations, we arrive at the final expression for the Casimir 
energy: 



^2*^ vac 4^^3 J 2n ^ 



^rcn 



5n 



UJ 



oo oo 

kdk ^ f , 

nr 



— 2 — tanh'^ xx 

A=±l 

= cosh 0+ tanh'^ >f+ + cosh 0- tanh'^ k 
= — 2iria, 



J cosh6'A , ,n ^ , AA' sinh^ 6'a' cosh e'_A' 

■( 1 ^ tt; — (1-tanh x-x) + > -, — , „ , 

I cosh6'_A ^ >irA'(cosh6'+ +cosh6i_) 

\ A' — i 1 

f cosh 0- tanh'^ k , 



(6.20) 

AA' sinh^ 6>a' cosh 6>_a' 1 , . 

rnsli , nnfih W ' ^ 



sinh6'± — , cosh0± = — 

K± K± 



(6 


21) 


(6 


22) 


(6 


23) 


(6 


24) 


(6 


25) 


(6 


26) 



where k, K±, lu, h± are now dimensionless, and all the square roots are taken in the algebraic sense, i.e., with a branch 
cut over negative real numbers and Rei/a; > Va; G C\(— oo,0). The integrand exponentially falls down with the 

increase of >f+ — \J + k"^, so the integral is convergent at infinity. It is easy to show that it is also convergent near 
zero. 

The integrand is complex-conjugate for the opposite values of then the result is real, as should be. Moreover, 
the only imaginary quantity in the integrand expression is 2ir]a^ so, to make the integral real, only even powers of rj 
should be present in the result. This means that the same integral as we have obtained above, applies to the case 
?7 < (although, in the above calculations, we assumed that r/ > 0). The rj parameter is present as a dimensionless 
combination rja. 

In the Maxwell {rj = 0) case, we obtain the following integral: 



vac 



1 



7) = 



kdk 



KldK+ 



(1 — tanh >ir)^ 
xtanh >c 



(6.27) 



where — y + k'^. Using the polar coordinates >(r,^, so that k = >ccos^, = >fsin^, ^ e [— 7r/2, 7r/2], and 
making the integration over ^, we obtain 



Tptcn 
jj2 vac 



1 



))=0 



Q , (1 — tanh >f)2 

K-^dx- — 

tanh K 



2 r(4)CR.(4) 
37r2a3 44 



5760a3 



720D^ ' 



(6.28) 



which leads to the classical result 



asimir 



r)=0 



9 1 pren 



ri=0 



2400-^ ■ 



(6.29) 



To find the corrections to this force, one should take partial derivatives of the integrand in the expression (|6.20p 
with respect to 77. This integrand, as seen, is symmetric as a function of K^,K_. Then one can easily demonstrate 
that 



r(=0 



dK. 



LU{K- 



(6.30) 



ri=0 



i.e., the partial derivative with respect to rj is proportional to the total derivative with respect to K^. The latter one 
vanishes after the integration over K^. The expression for the second derivative we are interested in is much more 
complicated, so we present it in a slightly different form, namely, when rj is set to zero after the differentiation, we 
use the polar coordinates x, ^ (see above) , since they are separated in the 77 = case. This second derivative is real 
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and has the following form: 
92 



a(?7a)2 



n=±i 



5n 



1 



77=0 



2>fsinh 2) 



{-2Ak^ coih2K + {^K- 7) cosh 4>< + 3 smh 4>tf + 16>^ + 7} - 8>tf + 4 
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2 cos 2^ f -2>!rsinh4>!r- cosh4>!r + 8>^2 -f 1 
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cos 4^ ■ 



1 - coth 2} 



{Ax{x+ l)smh4>^ + (4>^^ + Ax + i) cosh Ak ^ Ak^ - Ax-?,] 



2>tfsmh 2x 

After integrating over ^ e [— 7r/2, 7r/2], we obtain 

tt/2 



-■it/2 



n=±i 



ffn 



(6.31) 
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(144>!rsinh2>< + 7cosh6><- CdQx^ + 7) cosh2>!r) . 

15xsinh2 2>^^ ^ ' > 



Finally, the second derivative of the vacuum energy with respect to r\ reads 
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(6.32) 
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(6.33) 



Note that we have proved in Sec. IIVI that for |r/| < ir/Aa, the imaginary-energy solutions exist only for k = jr/j and 
they do not contribute to the vacuum energy which contains the integral over momentum k. Hence we can write E^™ 
instead of ReiJ^™ in the I77I <C 1/a case. In this case, taking (|6.28l) and (|6.33p together, we find the approximate 
expression for the Casimir energy and the Casimir force in the Maxwell-Chern-Simons electrodynamics: 
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(6.34) 
(6.35) 



which are precisely the same as those obtained using the zeta function approach (|5.27p . (|5.28p . However, here we have 
strictly accounted for the "quasi-zero" modes which were discussed qualitatively in the end of Sec. |Vl Our expression 
(|6.20p for the real part of the Casimir energy is exact for large \'q\D, i.e., at large distances. 



VII. DISCUSSION AND CONCLUSION 



Let us briefly discuss the results of our calculations. Originally, M. Frank and L Turan attempted to solve the 
problem we discussed here, using the Green's function method, in [33l |. But they have used a wrong identity, namely, 
(—92 -\-rf)e'^'^°'PdaAi3 — 0, in the MCS electrodynamics, that has reduced the dispersion relation for the MCS photon 
to that for a massive photon. The fact is, in 4 dimensions, the existence of the Chern-Simons term does not make 
the photons just massive, like it is in the 3-dimensional Maxwell-Chern-Simons electrodynamics j26| . instead, the 
photon possesses a more complicated dispersion relation (see Eq. p.lSp ). The result we have obtained differs both in 
sign and in magnitude from the one obtained in (33i] . It should be also mentioned that the calculation of the Green 
function within the MCS electrodynamics seems to be much more complicated than it was thought in [s^, so in the 
present paper we have used the two other methods, namely, the zeta function regularization and the summation and 
renormalization of the discrete sum involving the residue theorem. 

The zeta function analytical regularization is widely used for the calculations of the Casimir effect in various 
physical situations p^ . It automatically subtracts the vacuum energy of the infinite space (i.e., without the plates) 
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from the Casimir energy. In our calculations, we used the "true" zeta function regularization, in which the complex 
regularization parameter s, it depends on, controls the negative power of the frequencies in the series (jS.ip . Sometimes 
(see, e.g., [131 )i the space dimensionality d is chosen as the parameter of the analytical regularization, and, indeed, 
the expression for the vacuum energy depends on d through the Ricmann zeta function. In our case, it would not be 
strict enough, since the transformations we would make during the calculation of the vacuum energy (analogous to 
()5.17p . (IS.lSp . (I5.19|) ) would not converge together for any d E C The regularization with the parameter s avoids this 
problem, as it was mentioned in Sec. |Vl 

The method we have used to find the sum of a discrete series using the complex plane integral, is a kind of 
generalization of the Abel-Plana formula which is also widely used in Casimir effect calculations [41|. Indeed, one 
can generalize it to find the ex plic it integral expression for the series over the roots of a transcendental equation, 
including, e.g., Bessel functions [2^. The approach we have used in section IVll is another generalization of this type, 
however, it does not follow the approach developed in [25| . 

One should hold in mind that the calculations presented in our paper cannot provide a direct way to make experi- 
mental predictions. For the experimental purposes, we should take into account the finite conductivity of the plates 
and the dispersion of their conductivity, as well as some other aspects. Here we can only say how the presence of a 
small T] condensate, which violates Lorentz invariance, affects the dependency of the Casimir force on the distance 
between the plates. 

Again, we conclude that the correction is quadratic in rj and strengthens the attraction between the plates at 
relatively large distances of the order D < l/\ri\. The relative magnitude of this additional force compared to the 
Maxwell Casimir force increases quadratically with distance D, so large-separation Casimir effect measurements could 
give tighter constraints on 77. 

Recent observations [12] showed the 1 — 2% agreement between the experiment and the theory of the Casimir effect 
based on the conventional Standard Model, at distances of about R ^ 100. ..500 nanometer. Thus we can conclude 
that the leading ?7-correction to the Casimir force is less or about 1 — 2% of its value for rj = 0. This leads to the 
following constraint: 

< \/l% • (500nm)"i « 5 • lO^^eV. (7.1) 

Some experiments measure the Casimir effect at distances of several micrometers, with about a 10% accuracy (43}. 
and this can at least make the above constraint 3 — 4 times tighter. 

Though this constraint is very loose compared to those astrophysical and some other observations place on 77, it 
demonstrates a property of a quantum vacuum within the Extended Standard Model. Moreover, this is the first 
constraint placed on the Chern-Simons Lorentz-violating term based on the correct calculation of the Casimir effect 
in (3+1) dimensions. 
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